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Abstract We investigate strong-coupling properties of a two-dimensional ultracold Fermi 
gas in the normal state. Including pairing fluctuations within the framework of a T-matrix 
approximation, we calculate the distribution function n(Q) of Cooper pairs in terms of the 
center of mass momentum Q. In the strong-coupling regime, n(Q = 0) is shown to exhibit a 
remarkable increase with decreasing the temperature in the low temperature region, which 
agrees well with the recent experiment on a two-dimensional 6 Li Fermi gas [M. G. Ries, et. 
al., Phys. Rev. Lett. 114 , 230401 (2015)]. Our result indicates that the observed remarkable 
increase of the number of Cooper pairs with zero center of mass momentum can be explained 
without assuming the Berezinskii-Kosterlitz-Thouless (BKT) transition, when one properly 
includes pairing fluctuations that are enhanced by the low-dimensionality of the system. 
Since the BKT transition is a crucial topic in two-dimensional Fermi systems, our results 
would be useful for the study toward the realization of this quasi-long-range order in an 
ultracold Fermi gas. 

PACS numbers: 03.75.Hh, 05.30.Fk, 67.85.Lm. 


1 Introduction 

The advantage of an ultracold Fermi gas is the high tunability of various physical param- 
etersia. Using a tunable pairing interaction associated with a Feshbach resonance, we can 
now study a Fermi superfluid from the weak-coupling regime to the strong-coupling limit in 
a systematic manner—. In addition, the system dimensionality can be also tunable by using an 
optical lattice potential™. Indeed, a two-dimensional Fermi gas has recently been realized by 
using this techniqu e 4 i 5 i 6 i 7 i 8 i 9 . Because of these experimental developments, strong-coupling 
properties of a two-dimensional Fermi gas has become an interesting and realistic research 
topic in cold Fermi gas physic s t 0 ] 1 1,i2 1 i3 1 i4 1 i5 i 1 6i 1 7,i^ 
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In contrast to a three-dimensional Fermi gas, the superfluid long-range order is prohib¬ 
ited in the two-dimensional case, because it is completely destroyed by low dimensional 
superfluid fluctuations' —< 2 2 . However, a two-dimensional Fermi gas is known to be able to 
still exhibit superfluid properties, when the Berezinskii-Kosterlitz-Thouless (BKT) phase 
transition occur s 2 *i 22 . 

Recently, the observation of this quasi-long-range order was reported in a two-dimensional 
6 Li Fermi gas—. While the BKT transition is theoretically explained on the viewpoint of 
vortex-antivortex pair annihilation, this experiment determines the BKT phase transition 
temperature T^ r as the temperature below which the number Nc(Q = 0) of Cooper pairs 
with zero center of mass momentum (Q = 0) remarkably increases— Since the vortex- 
antivortex pair annihilation is not observed in this experiment—, it is a crucial issue to check 
whether or not determined from the temperature dependence of Nq(Q = 0) can unam¬ 
biguously be identified as the BKT transition temperature predicted theoretically. 

In this paper, to examine this, we investigate a two-dimensional Fermi gas near 7^/. We 
discuss how to evaluate Nc(Q = 0) in a strong-coupling T-matrix approximation (TMA). 
Although this strong-coupling theory cannot describe the BKT transitio n .— ^ , we show that 
the observed remarkable increase of this quantity below a certain temperature (which is 
experimentally identified as the BKT transition temperature 8 ") can be explained without as¬ 
suming the BKT phase transition. We also present an alternative explanation for this phe¬ 
nomenon on the viewpoint of strong pairing fluctuations that are enhanced by the low- 
dimensionality of the system. Throughout this paper, we take Ti = = 1 and the two- 

dimensional system area is taken to be unity, for simplicity. 


2 Formulation 


We consider a two-dimensional uniform Fermi gas consisting of two atomic hyperfine states, 
described by the BCS Hamiltonian, 

H = 52 ^>P c p,a c p,o ~ U C p+q/2,-\ C -p+q/2,l C -p'+q/2,l c p'+q/2 ,f- ( 1 ) 

P ° p p' q 


Here, c pa is a creation operator of a Fermi atom with pseudospin a = j . J, and two-dimensional 
momentum p = ( p x ,Py )• %p = P 2 /(2m) — p is the kinetic energy, measured from the Fermi 
chemical potential p, where m is an atomic mass. The pairing interaction —U (< 0) is as¬ 
sumed to be tunable by adjusting the threshold energy of a Feshbach resonance. As usual, we 
measure the interaction strength in terms of the two-dimensional i-wave scattering length 
fl2D> which is related to U asa£ 1/1/ = (m/ 2n) In (^f^d) + Y, P >k F m/p 2 (where kp = \f2nN 
is the Fermi momentum, with N being the total number of Fermi atoms). 

Many-body corrections to Fermi single-particle excitations can be conveniently de¬ 
scribed by the self-energy E(p,ico n ) in the single-particle thermal Green’s function, 


G(p,i(Pn) = T 


To,.. Cp E(pJco n ) 


(2) 


where ft>„ is the fermion Matsubara frequency. In the T-matrix approximation, the self¬ 
energy E(p,ico n ) is diagrammatically described as Fig. [3 which give s— 


E(p, i(O n ) = T £ r(q,iv n )Go(q — p, iv n — i(O n ). 

qjVn 


( 3 ) 
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Fig. 1 Self-energy Z(p,ia>„) in the T -matrix approximation (TMA). The solid line represents the non¬ 
interacting Green’s function Go, and the dashed line denotes the pairing interaction —U. 


Here, v„ is the boson Matsubara frequency, and Gq 1 (p, ico„) = ico n — E, p is the single-particle 
Green’s function in a free Fermi gas. The particle-particle scattering matrix r(q,iv n ) in 
TMA has the for m—a 7 . , 

r(q,iv n ) = - l _ fj V / _ x , (4) 


Un(q.iv n ) ’ 


where 


n{q,iv n ) = T £ Go(p + | ; *v„-Hffl„)G 0 (-p+|,-iffl„) 


(5) 


P,1(0,1 


is the lowest-order pair-correlation function, describing fluctuations in the Cooper channel. 
The equation for the total number N of Fermi atoms in TMA is given by, 


N = 2T £ G{p,i(O n ). 

p,ia>„ 


(6) 


This number equation may be divided into the sum of the free-fermion contribution, 

Nq = 2T £ Go(p,ia>„), (7) 


p,l(On 


and the fluctuation correction 8N described by the TMA self-energy E(p, i( 0 n ) in Eq. (O, 


8N = 2T £ [G(p, icon) -Go(p, i(On)]. 

p.iio,, 


( 8 ) 


Pairing fluctuations in the BCS-BEC crossover region are physically understood as the rep¬ 
etition of the formation and dissociation of preformed Cooper pairs, that eventually become 
tightly bound molecular bosons in the strong-coupling BEC limit. Keeping this in mind, and 
writing Eq. © as 5N = 2 Y,Q n (Q)i one may regard 

n(g) = r£r(g,/v„)T £ G 0 {Q^p,iv n -i(O n )Go{p,i(O n )G(p,i(On) (9) 

iV n P,1(0,, 

as the number of preformed Cooper pairs with the center of mass momentum Q. Indeed, 
in the strong-coupling limit (where p —> — °°), n(Q) is reduced to the ordinary momentum 
distribution in an ideal Bose gas, 


Q 2 


2 M 


— Mb , 


n(Q) = n B 

where «b(®) is the Bose distribution function, M = 2m is a molecular mass, and 

J UB=2Mln(-=^) (<0) 

V ^bind / 


( 10 ) 


(ii) 
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Fig. 2 (Color online) (a) Calculated momentum distribution function n(Q) of Cooper pairs with respect to 
the center of mass momentum Q. We take ln(A:F^ 2 D) — —0.59. The dotted line shows the approximate result 
in Eq. {I3j. 


is interpreted as the Bose chemical potential, with Z?bind = l/(ma? D ) being the binding en¬ 
ergy of a two-body bound state. Although such a molecular picture gradually becomes worse 
as one approaches the weak-coupling BCS regime, n{Q) in Eq. Q is still a useful quantity 
to grasp the bosonic character of the system in the BCS-BEC crossover region. In this pa¬ 
per, thus, we identify this momentum distribution function n(Q ) with the observed number 
Nq ( Q) of Cooper-pair bosons with the center of mass momentum Q in a two-dimensional 
6 Li Fermi gas&— 

We briefly note that the two-dimensional T -matrix approximation we are using in this 
paper does not give the BCS-type superfluid phase transitio n 1 ? . '— , and also cannot describe 
the BKT phase transitio n 14 i 23 ' 24 . Thus, the comparison of our TMA result on n(Q) with 
the recent experiment—^ provides a useful information about whether or not the observed 
anomalous increase of Nq(Q = 0) is really a clear signature of the BKT transition. 

In this paper, to compare our TMA results with the recent experiment on a two-dimensional 
6 Li Fermi gas^a&, we take Inikpa^o) = —0.59. At this interaction strength, we first calcu¬ 
late the chemical potential p(T) from the number equation <(6]), and then evaluate n(Q) in 
Eq. <[9} at various temperatures. We briefly note that Ref.— reports the BKT phase transition 
temperature Tg^. = 0.1297p at this interaction strength, where 7p is the Fermi temperature. 


3 Momentum distribution function of Cooper pairs 


Figure [2] shows the momentum distribution function n(Q) in the low temperature region 
where the BKT phase transition was experimentally reported (T^j/Tf = 0.129). In this fig¬ 
ure, n{Q) in the low momentum region is found to be remarkably enhanced with decreasing 
the temperature. To see to what extent this behavior reflects bosonic character of the sys¬ 
tem, it is convenient to approximately evaluate the TMA self-energy E(p, ico n ) in Eq. ([5} by 
employing the so-called static approximation as- 27 ’22, 


E(p. ia>„) 


T £ r(q,iv n ) 

q,iv n 


Go{-p..—i(O n ) 


-Ap G G 0 (-p,-i(O„). 


( 12 ) 


Here, Ape is the so-called pseudogap parameter in the literature)— This approximation as¬ 
sumes that pairing fluctuations described by r(q,iv„) are enhanced in the low-momentum 
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Fig. 3 (Color online) Calculated Fermi chemical potential as a function of temperature, when In (&f 02 d) = 
—0.59. ftTh is the chemical potential which satisfies the Thouless criterion—. Although U is very close to /ixh 
when T/Tp < 0.15, the former is always slightly lower than the latter. 


and low-energy region. In addition, as shown in Fig. [3] the Fermi chemical potential jl is 
negative and |mI/ £ f > 3 when ln(fcpa 2 D) = —0.59 (where £p is the Fermi energy), indi¬ 
cating that the system at this interaction strength is already in the strong-coupling regime. 
Including this, one can approximate n(Q) in Eq. |9} to 


_2£f^” B (y2M MB )’ (13) 
^bind 

where we have used the strong-coupling expression for the pseudogap parameter, Ape — 
2 a/Ef [£f + -Ebind /2] , as well as fi ~ — Fbind/2, in obtaining the last expression (where the 
two-body binding energy Ebind is given below Eq. 11 111). As shown in Fig.(2j Eq. 1131 ) well 
describes the TMA momentum distribution function n{Q ) when T/Tp =0.11. This clearly 
indicates that the behavior of n(Q) shown in Fig.Qstrongly reflects the bosonic character of 
this system^i, originating from the formation of preformed Cooper pairs. 

Figure [4] shows n(Q = 0) when ln(kpa 2 p>) = —0.59. In this figure, we find that the 
remarkable increase of this quantity around Q = 0 seen in Fig. [2] starts to occur when 
T — T^ t = 0.1297p. We also find that this temperature dependence agrees well with the 
recent experiment on a two-dimensional 6 Li Fermi gas>£. As mentioned previously, since our 
TMA gives no superfluid phase transition around 7^^ = 0.1297 f, this agreement indicates 
that the observed remarkable increase of the number of Cooper pairs with zero center of 
mass momentum does not necessarily mean that the system is in the BKT phase. Further 
experimental studies would be necessary to confirm that the BKT phase is really realized 
below = 0.1297p at this interaction strength. 

The enhancement of n{Q = 0) in TMA is due to strong pairing fluctuations enhanced 
by the two-dimensionality of the system. To explain this, we expand the denominator of the 
particle-particle scattering matrix F (q, iv n = 0) in Eq. 0 around q = 0 as 


"(G) = 


— 2M(M + y / M 2 +^p G ) 


A 2 

^PG 


n B 


2^ ~ Mb ) = 


1 + 


r ( q .,v„ 0) ~ [l _ U n(q = 0,iv n = 0)]+aq 2 ’ 


(14) 


where a is a positive constant. When the so-called Thouless criterion^, 1 — [7/7(0,0) = 0, 
is satisfied (which is known as the condition for the BCS-type superfluid phase transition 
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Fig. 4 (Color online) Calculated momentum distribution function n(Q) at Q = 0 when ln(A:p^ 2 D) — —0.59. 
We also show the recent experiment on the number of Cooper pairs with zero center of mass momentum in a 
two-dimensional 6 Li Fermi gas£. The two dashed lines are experimental fitting lines. Reference! determines 
the BKT phase transition temperature = 0.1297p as the temperature at which these two lines crosses 
each other. Each of theoretical and experimental result is normalized by the value at T^ t . 


in the three-dimensional case), the self-energy in Eq. dT2l) diverges in the two-dimensional 
case, reflecting the logarithmic divergence of the momentum integration of Eq. ® [~ 
fqdq(l/q 2 )] 23 i 24 . As a result, the chemical potential jx never reaches the value (= jUxh) 
that satisfies the Thouless criterion. (See Fig. [3]) Because of this, TMA does not give the 
superfluid phase transition in a two-dimensional Fermi gas. Flowever, in the low temperature 
region where n ~ /Xjh in Fig. [3] one expects that 1 — UFI (0,0) ~ 0. In this case, since the 
particle-particle scattering matrix in Eq. 0 also has the meaning of an interaction between 
Fermi atoms, this effective pairing interaction is found to remarkably be enhanced by many- 
body effects in the low-energy and low-momentum region, which positively contributes to 
the formation of Cooper pairs. As a result, n[Q = 0) in TMA anomalously increases in this 
temperature region, as shown in Fig.|4] 

Before ending this section, we briefly note that Ref. 8 reports that Tg^. is higher for 
a weaker pairing interaction, which is opposite to the theoretical prediction that the BKT 
phase transition temperature decreases with decreasing the interaction strength I S*™ . In this 
regard, the temperature dependence of n(Q = 0) in TMA exhibits the same tendency as 
the experimental result— (although we do not show the result here), which also implies the 
importance of pairing fluctuations in considering this physical quantity. 


4 Summary 

To summarize, we have investigated a two-dimensional ultracold Fermi gas in the BCS-BEC 
crossover region. Including pairing fluctuations within a T -matrix approximation (TMA), 
we calculated the momentum distribution function n(Q) of Cooper pairs in the normal state. 
We showed that the anomalous enhancement of this quantity at Q = 0 that have recently 
been observed in a two-dimensional 6 Li Fermi gas£ may be quantitatively explained with¬ 
out assuming the BKT phase transition. That is, even in the normal state, strong pairing 
fluctuations in the BCS-BEC crossover region that are enhanced by the two dimensional¬ 
ity of the system also lead to the enhancement of n(Q = 0) around the temperature which 
was experimentally identified as the BKT phase transition temperature T^ r . This indicates 





Title Suppressed Due to Excessive Length 


7 


that further experimental studies are needed to confirm that the system is really in the BKT 
phase below 7^^,. Since the BKT transition is one of the most crucial topics in the field of 
two-dimensional Fermi superfluid, our results would contribute to the study of this quasi- 
long-range order in cold Fermi gas physics. 
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